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A self consistent solution to Dirac equation in a Kerr Newman space-time with M 2 > a 2 + Q 2 is 
presented for the case when the Dirac particle is the source of the curvature and the electromagnetic 
field. The solution is localised, continuous everywhere and valid only for a special choice of the 
parameters appearing in the Dirac equation. 
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. Bound solutions of the Einstein Maxwell Dirac equations are of interest as a link between quantum mechanics and 
general relativity in a regime where both theories are well tested. Finster, Smoller and Yau ply, Kg along with other 
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collaborators showed that no time periodic solutions exist for the Dirac equation in Schwarzschild, Reissner Nordstrom 
and Kerr Newman (M 2 > a 2 + Q 2 ) space-times. Dafermos |4J has confirmed this result by another method. On the 
other hand Schmid Q showed that such time periodic solutions to the Dirac equation do exist for the extreme Kerr 
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In these investigations, it was assumed that the Dirac particle mass m e is very small as compared to the black 
CNJ . hole mass M, that is, it is an external test particle. I show in the self consistent case, when the Dirac particle is the 
source of the curvature and the electromagnetic field, that a solution is possible for certain values of the parameters 
(m e , Q, a) appearing in the Dirac equation. Further this solution is continuous everywhere and normalisable. 

Such a solution corresponds to a generalisation of the free particle Dirac equation in Minkowski space to include 
the effects of the curvature produced by the particle. The ordinary free particle solutions of the Dirac equation 
are completely delocalised; the curvature however, now causes the Dirac wave functions to be localised over a region 
£f~) \ comparable in dimension to the Compton wavelength of the particle. Note that the wave function still has an enormous 

■ spread in comparison to the dimensions of the event or Cauchy horizons of the particle. 
""q I As is usually the case in single particle Dirac theory we ignore the self interaction between the particle and its 
own electromagnetic field. Further we assume that we need only to look for solutions with zero azimuthal angular 
momentum about the rotation axis. In the next section the equations resulting from Chandrasekhar's separation of 
the Dirac equation in a Kerr Newman spacetime 0, are simplified by applying these assumptions. 

A necessary requirement for a wave function is that it be single valued at a point. We carry out a bifurcation 
analysis of the Dirac equations for 9 and r to find their bifurcation points. The requirement that these lie outside 
the physical region of interest in both the coordinates 9 and r gives us conditions on the parameters appearing in the 
Dirac equation. Here we confine ourselves to the simplest case. The equations are then solved and it is found that 
the angular part of the solution is continuous and normalisable. The radial equation is also solved in each of the three 
regions (i) the exterior, (ii) the ergosphere between the event and the Cauchy horizons and (iii) the region inside the 
Cauchy horizon containing the singularity. The resulting radial solutions are continuous in all three regions and so 
a simple matching of the wave functions at the horizons gives us a unique normalisable solution, with zero energy 
(u) = 0, A = 0). 

In the last section, the nature of this solution is discussed. In particular, if the fields are centered at r — 0, it is not 
obvious what it means to have a finite probability of finding the source at r ^ 0. I suggest that this can be explained 
in a manner similar to the explanation for the Klein paradox of Dirac theory. Namely, in the region of the order of the 
Compton wavelength within which the wave function peaks, the field strengths are so large that virtual pair creation 
cannot be neglected. The antiparticle thus produced can recombine with the "original" particle at r — leaving us 
to detect the particle now at r ^ 0. Consider next the case of the only allowed energy eigenvalue zero. For a single 
bound particle in a potential, the energy can of course be negative or zero. 
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II. SOLUTION OF THE DIRAC EQUATION 

Following the notation of Chandrasekhar's book 7], the Kerr Newman line element in Boyer Lindquist coordinates 
is given by 

(ds) 2 = ^(dt-asm 2 9dcj)) 2 -J:(^-+d9 2 )-^-{adt-(r 2 +a 2 )d^) 2 , (la) 

S = r 2 + a 2 , (lb) 
A = r 2 -2Mr + a 2 + Q 2 . (lc) 

As we are using natural units, the mass appearing in the metric (M — > ^4f-) has a size comparable to the horizon 
length. Whereas, the mass m e which appears in the Dirac equation is (m e — > 2 ¥ s ) has a size which is the inverse of the 
Compton wavelength. Thus we use different symbols to distinguish the two length scales in the problem. They both 
of course refer to the same mass but differ according to which of the universal constants G, fi, c was being suppressed. 

The Dirac equation can then be separated by writing the four components of the spinor wave function in the form 

tfi - 2 - a-e*^*-"*), (2a) 

r — ia cos 9 

* 2 = R$Si e ,(m ^ ut) , (2b) 

*3 - R k S -h e i(m0 ~ wt) , (2c) 

* 4 = J-J-e'^""''. (2d) 

r + ia cos 9 

I have slightly modified the notation to use u> for the frequency rather than er. The R±± are functions of r only 
while S±i are functions of 9 only. If this form is substituted into the Dirac equations, we get the following angular 
equations. 

— S^i = — (mcsc#H — cot 6)S, i — (A — am e cos 9)S_i, (3a) 
up 2 2 2 2 

-jtS_i — (mcsc9 — i cot9)S_i + (A + am e cos0)5 + i . (3b) 
ap 2 2 2 2 

(3c) 

We now set m = and carry out a bifurcation analysis of these equations. The angle equations © will simplify 
under the substitutions 



r ±4 = Vsin0S ± i, (4) 



to yield 



-^T+i = - (A -am e cos 6)T_i, (5a) 
ac' 2 2 

d 

(5c) 



jT_i = (A + am e cos 6»)T +i . (5b) 



These have bifurcations at 



(A) 2 = (am e cos 6») 2 . (6) 



Our condition that the wave functions be single valued means that we wish to have no bifurcations in the range 
< 9 < 7T. This is possible if and only if |A| > am e or A = 0. The second case corresponds to an s wave orbital state 
in the usual separation for the Dirac equation 9] in a central potential. I have only considered this case here. The 
angle equations then are 

—T,i = am P cos9T i, (7a) 

d0 +2 -5' v I 

—T i = am e cos 6T,i. (7b) 

d9 ~2 +2 v > 
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These are trivial to solve and we obtain 



Si = A(e am " sine ± e - amcSine )/(V2iin0), (8a) 



S_i = A(e am ^ sine T e- am ' siDe )/(V2^9). (8b) 

These solutions while they have a formal singularity at the poles give a finite probability density every where. 

Next consider the radial equations. Again I find it convenient to use a slightly different variable from these used 
by Chandrasekhar (0), we define R i = y/ARi 

d i , , uj(r 2 + a 2 ) + Xa / , . A — im e r „ . , . . 

= — ^ + — ^fi-iW- ( 9a ) 



d . , . u(r 2 + a 2 ) + Aa A + im e r > 

—R_i(r) = —i t R_i(r)-\ = — R .i(r). (9b) 

dr = w A s w a/A ' 



The corresponding bifurcation points lie at those values of r which satisfy 

{iu{r 2 + a 2 ) + \a) 2 = (A 2 + m e r 2 )A. (10) 

This fourth order equation permits a wide range in the 5 dimensional parameter space (m e , A, Q, lo, m) of the Dirac 
equation for which the bifurcations lie outside the region of interest. It is worth remarking that such a solution must 
always have uj 2 > m 2 . Further the solutions with lo = ±m e are unconditionally stable. As we are considering only the 
case of no orbital angular momentum, A = 0, equation l|l(J|l also permits solutions with lo = which have bifurcations 
at r = M ± VA and r = 0. I will consider only this last case (to = 0) in this work. For lo = and A = 0, the radial 
equations (JSJ simplify to yield 

= -^-R-i(r), (Ha) 



r = -7^fl + i(r). lib 
dr 2W VA 2 



Again as in the case of the angle equations, these can be solved quite easily. Consider first the exterior region 
r > M + \/~K if we change the independent variable to u = \/A + M In (r — M + VA) the radial equations 111(1 
simplify to yield 



d 

du 
d 



—R + i(u) = —im e R_i{u) 1 (12a) 



— R_i(u) = im e R,i(u). (12b) 
du 2 T2 



We can thus write down the solution in the exterior region which is finite as r — > oo as 



g— TO e VA 

%/Ai?! = B -= -, (13a) 

[r-Af + VA] m = M 

„— m e \/A 

= -tS == . (13b) 



Next consider the ergospheric region between the event horizon r+ = M + \J~K and the Cauchy Horizon r_ = M—^/A. 
Here we transform the independent variable to v = A| + Mln (r — M + an d obtain the solution 

y/ARi = C*e™ c Vl A l[ r _ M + iJ\A\] m ° M + D = , (14a) 

2 V [ r _ M + i y /m]m c M' 

R i = iCe mi 'V\ A \\r - M + iy / \A\] m " M - iD = . (14b) 

" 5 [r- M + 2 y|A|] m " M 
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Finally in the region interior to the Cauchy horizon r < r_ we use the new independent variable w = vA + 
M In (r - M + VA) to obtain the solution 

VARi = Ee m ^[r - M + \/A] m « M + F == , (15a) 

p -m e \fA 

R i = -iEe m ^\r-M + VA] m ^ M + iF = . (15b) 

Thus the solutions are continuous in all three regions, the only effect of the bifurcation at r+ and r_ being an ambiguity 
as to whether it is my choice in ea. l|l 41 which is correct or its complex conjugate which is to be used in this region. 
As the solution is continuous in all three regions, the coefficients B,C, D, E, F can be found up to a normalisation 
constant by requiring that the solutions be continuous at the horizons. So we finally obtain the solution 

ARi = B = , (r > r+) (16a) 

" [r-M + VA} m * M ~ 

e -im e y/\A\ 

ARi = B = , (r_ < r < r+) (16b) 

5 [r- M + is/\A\] m " M " " 

VARi = Be m ^[r - M + y/A} m ? M , (r < r_) (16c) 

R i = -iB = , (r>r+) (16d) 

" 3 \r-M + VN m ° M ~ 



/\A\ 

R_i = -iB (r_ < r < r + ) (16e) 



[r - M + iy f \K\] m ' M 
R_i = -iBe m ^[r-M + y/A] m ' M (r < r_). (16f) 



III. DISCUSSION 



The solution we have just obtained is normalisable and thus represents a proper Dirac spinor wave function. There 
is of course a second possibility due to the bifurcations at the horizon, which is to take the complex conjugate form 
inside the ergosphere. This however will give rise to the same Dirac currents etc., outside and is thus indistinguishable 
form the form presented here. 

Next consider the question of the only allowed energy eigenvalue, zero. In any bound state problem, we expect 
the energy of the state to be less than or equal to the asymptotic value of the binding potential. In this case if we 
consider the curvature to be producing the attractive force localising the particle close to r ~ 0; then the zero energy 
eigenvalue should evoke no surprise. 

Finally there is the question of what is meant by finding the source of the field at a point other than where the 
field originates from. Again if this is considered as a potential problem for the Dirac field this can be easily explained. 
The wave function given by eq. l|16|) is localised over a region of the order of the Compton wavelength of the particle. 
In this region the field strengths are very high and so virtual pairs of the same particle type have finite amplitude of 
appearing. The antiparticle of such a pair can then combine with the " original" particle leaving the now real particle 
of the pair to be found in the position r =/= 0. 

While the present solution might appear overly simplistic, it is the fact that it exists which is important. This 
suggests that other such solutions may also be found in the rest of the allowed parameter space which was not explored 
here. 
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